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Abstract 

Short-period binaries represent extreme cases in the 
generation of stellar coronae via a rotational dynamo. Such 
stars are important for probing the origin and nature of 
coronae in the regimes of rapid rotation and activity satu- 
ration. VW Cep (P = 0.28 d) is relatively bright, partially 
eclipsing, and very active object. Light curves made from 
Chandra/HETGS data show flaring and rotational modu- 
lation, but no strong eclipses. Emission lines are broader 
than instrumental, indicating emission from both binary 
components. Velocity modulation of emission lines is in- 
dicative of geometric structure of the emitting plasma. 
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spheric spots and coronal emission) ; Doppler image maps 
of VW Cep ( |Hendry fc Mochnacki 2 000 ) showed large po- 
lar spots. 

IChoi fc Dotani 1998l have analyzed ASCA spectra, and 
found a flux of about 1 x 10~ n ergscm _2 s _1 , and two 
component model temperatures of 7 and 22 x 10 6 K (log T ~ 
6.8 and 7.3) with about equal emission measures. They 
obtained significantly reduced abundances of Fe, Si, Mg, 
and O, but a Solar value for Ne. A flare also occurred 
during the observation, with a factor of three increase 
in the count rate. They used the flare emission measure 
and loop-scaling models to estimate a density of about 
5 x 10 10 cm~ 3 . The Ginga observations IjTsuru et al. 1 992) 
showed a thermal plasma temperature in excess of 10 s K, 
a flux similar to that determined bv lChoi fc DotanT~1 998 
no rotational modulation, and Fe K flux lower than ex- 
pected. 



1. Introduction 

VW Cep (HD 197433) is a W-type W UMa binary — a 
contact binary in which the more massive and larger star 
has lower mean surface brightness. The primary (deeper) 
photometric eclipse is the occultation of the smaller star. 
VW Cep has an 0.28 d (24 ks) orbital period. It is par- 
tially eclipsing (i = 63°), with component spectral types 
ofKO Vand G5 V l|Hill 1989l|Hendry fc Mochnacki 2000| |. 
It is one of the X-ray-brightest of contact binaries. 

Among the coronally active binaries, activity is strongly 
correlated with rotation l|Pallavicini 1989| . but saturates 
at periods below one day IjVilhu fc Rucinski 1983l|Cruddace 
The trend has been referred to as "super-saturation" since 
instead of reaching a plateau, the activity level actually de- 
creases with increasing rotation rate ( Randich 1998 , J ardine 
IJames et al. 2000irSt§pien, Schmitt fc Voges 2001| ). IJardine 
argued that super-saturation occurs because loops are large 
and are unstable to the coriolis forces as they exceed the 
co-rotation radius: extended loops are swept to the poles. 
On the other hand, |St§pien, Schmitt fc Voges 2001| con- 
tended that loops are quite compact (relative to the stellar 
radius), but that the dynamics of surface flows in contact 
systems conspire to clear equatorial regions. The two sce- 
narios are similar, but they differ in an important respect: 
X-ray sources are either low volume and dense, or large 
volume and rarefied. Optical light curve modeling is con- 
sistent with either (assuming correlation between photo- 



2. High Resolution X-Ray Spectra 

We observed VW Cep for 116 ks with the Chandra High 
Energy Transmission Grating Spectrometer (HETGS) on 
August 29-30, 2003 (observation identifier 3766). The in- 
strument has a resolving power (E/AE) of up to 1000 
and wavelength coverage from about 1.5 A to 26 A in two 
independent channels, the High Energy Grating (HEG), 
and Medium Energy Grating (MEG). 

VW Cep has also been observed in X-rays with the 
Chandra Low Energy Transmission Grating Spectrometer 

^n"?""™^ :: ^ e ^ ^ >u P ree 2 003|) , and with the XMM- 

rlyfbservatorv HGondoin 20040 . The only other 

contact binary observed with HETGS is 44 Boo ( |Brickhouse, Dupr 



Ohsnrvatinns with the HETGS provide spectral resolution 
,V Lin uli -J.J.U1 re q U j rec j t better determine line fluxes, 



*w aVt?feiLg Liis^ ^Lnd profiles. Ultimately, we wish to deter- 
mine which coronal characteristics are truly dependent 
upon fundamental stellar parameters. Combination of VW Cep 
data with that from LETGS and XMM-Newton, as well as 
with results for other short period systems (e.g., 44 Boo, 
ER Vul) will help us to understand the basic emission 
mechanisms. 

Here we concentrate on preliminary HETGS results. 
Data were calibrated with standard CIAO pipeline and 
response tools 1 . Further measurements and analysis were 



1 |http : // cxc . h arvard . edu/ c i ao 



Proc. 13th Cool Stars Workshop, Hamburg, 5-9 July 2004, F. Favata et al. eds. 
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Figure 1. HETGS spectrum of VW Cep, 116 ks exposure. MEG (16433 counts) is the black line, and REG (5269 
counts) the gray. Insets show detail of the REG spectrum (left) and MEG O VII triplet and N VII region (right). 



done in ISIS (Houck & Dcnicola 2000) and with custom 
code using ISIS as a development platform. 

3. The Spectrum 

Figure ^ shows the counts spectrum obtained from the 
HETGS full exposure. The spectrum is qualitatively typ- 
ical of coronal sources: a variety of emission lines from 
highly-ionized elements formed over a broad temperature 
region, from Ovu, Nvn, Neix, and Fexvn (logT ~ 6.3- 
6.7), up to high temperature species like S XV, S XVI, Caxix 
(perhaps), and Fexxv (logT ~ 7.2-7.8). It is also appar- 
ent that iron has a fairly low abundance relative to neon, 
given the relative weakness of the 15A and 17A FeXVH 
lines relative to the Neix 13 A lines. The observed flux in 
the 2-20 A is 7 x 10" 12 ergs cm s , and the luminosity 
(for a distance of 27.65 pc) is 6.5 x 10 29 ergss _1 . 

Cursory inspection of the density-sensitive helium-like 
triplet lines (Ovu, Neix, Mgxi) does not reveal ratios 
obviously far from their respective low-density limits. 

4. Emission Measure 

The differential emission measure (DEM) is a one dimen- 
sional characterization of a plasma, and can be defined as 
iV 2 dV/dT, in which N e is the electron density, V is the 
emitting volume, and T the temperature. The DEM is an 
important quantity because it represents the radiative loss 
portion of the underlying heating mechanism. An emission 
measure can be derived from measurements of line fluxes 
and some assumptions about the homogeneity and ion- 
ization balance of the emitting plasma. Derivation of the 
emission measure thus relies on detailed knowledge of fun- 



damental atomic parameters. Even given perfect knowl- 
edge of the ionic emissivities, their contribution functions 
versus temperature are broad so the emission integral can- 
not be formally inverted. Hence there are many meth- 
ods to regularize the problem and solve for the emission 
measure and elemental abundances. For a first estimate, 
we use the simple method described bv lPottasch 19631 in 
which the DEM is approximated by a ratio of line lu- 
minosity to average line emissivity at the temperature of 
the maximum emissivity. Since there is a degeneracy in 
elemental abundance and DEM, the relative abundance 
of one element is scaled to bring their DEM into bet- 
ter agreement with the locus of some other element with 
peak emissivities near the same temperatures (typically 
Fe, since it has many ions over a broad range in temper- 
ature). 

Figure [5] shows our provisional DEM, using the As- 
trophysical Plasma Emission Database (APED) for line 
emissivities (Smi trTet al. 2001^1 . the ionization balance of 
Mazz otta et al. 19981 and Solar abundances o flAnders fc Grevesse 
There is large scatter in the values; the solid curve shows a 
smoothed average. There is sharp structure, with a fairly 
large peak at logT = 6.9. This DEM can be used to 
generate a synthetic spectrum for further refinement of a 
plasma model. We found that with this distribution, we re- 
quired a neon abundance of 1.0 (relative to Solar), and 0.3 
for iron, confirming a trend already inferred from visual 
inspection of the spectrum. In iterative DEM reconstruc- 
tion techniques, sharper structure can be modeled than 
in the simple method. The fact that the simple method 
shows a strong peak hints at a highly structured DEM in 
VW Cep. 
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Figure 2. A provisional emission measure constructed from 
the ratio of line fluxes to mean emissivity over 0.3 dex. 
Points are plotted for each line used, and different ele- 
ments points were scaled as a group to be roughly consis- 
tent with others. The solid line is a smoothed average of 
the individual points. The integrated emission measure is 
about 6.8 x 10 52 cm~ 3 . 

5. Light and Phase Curves 

Given the short orbital period, the observation covered 
almost five revolutions without interruption. Such phase 
redundancy is important to discriminate intrinsic variabil- 
ity from that caused by rotational or eclipse modulation. 
Figure |3] shows light curves in several wavelength regions 
and a hardness ratio. There was much variability. The 
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Figure 3. Count rates of VW Cep in 2 ks bins. In the upper 
panel are light curves extracted from all diffracted photons 
(1.5-26 A; upper heavy (black) curve), a "soft" band (12- 
26A; lower, thinner curve (red) ), and a "hard" band (1.5- 
8.3A; lower, thicker curve (blue)). The lower panel shows 
a hardness ratio (solid histogram), and the median of the 
hardness for the first 70 ks (light dotted line). 

hardness shows that the large increase near after 80 ks is 
a flare, by definition of being hotter: proportionally more 
flux was emitted at shorter wavelengths, which are very 
sensitive to high temperatures via the thermal continuum 




Phase 



Figure 4. The phase-folded X-ray count-rate curve of 
VW Cep, excluding flare times. Bins are A(j) = 0.025, 
or 600 s, with a cumulative exposure per bin of 1800 to 
2400 s, depending on the phase. 

emission. Conversely, the bump in count rate between 20- 
30 ks does not show in hardness, and must be due to 
rotational modulation. 

If we exclude the flare times and phase-fold the re- 
maining data (about 3 orbital periods), we easily see the 
rotational modulation. Figure 01 shows this curve, using 
the ephemeris of Pri bulla et al. 20021 No eclipses are ap- 
parent, but instead, quasi-sinusoidal modulation of about 
18%. The visual light curve ( |Pribulla, Parimucha fc Vanko 2000| ) 
shows a similar amplitude (16%), but is very different 
qualitatively, with strong minima at phases 0.0 and 0.5, 
and continuously variable in between (a trademark of W UMa 
systems) . Without additional information, the X-ray light 
curve is not sufficient for localizing emission to one star or 
another. All we can say is that there is some asymmetric 
distribution, and possible occultation of large structures. 

6. Velocity Modulation 

In the high-resolution spectrum, additional information 
is available in the line shapes and centroids. The lines 
of VW Cep are broader than the instrumental resolution 
when integrated over the entire observation. This means 
that we are most likely sensitive to orbital Doppler shifts 
of one or both components, or to rotational broadening. 

In order to investigate further, we have measured the 
mean wavelength of events in the best exposed line, Nex 
(12. lA), in phase bins of 0.1 (also excluding flare times). 
The one-sigma envelope is shown in Figure It roughly 
follows the radial velocity of the primary (more massive) 
star (heavy dark sinusoid), except near phase 0.5, where 
it does a flip-flop, moving red-ward, then blue-ward. Our 
current hypothesis is that we are seeing a bias in the 
centroid during the transit of the secondary which first 
blocks the approaching limb of the primary, then the re- 
ceding limb as the transit progresses. The dashed sinu- 
soidal curves in the figure denote approximate photosphcric 
rotational velocities on opposite sides of the primary. The 
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Figure 5. Mean wavelengths of photons in the Ne x 
(12.1 A) region of the spectrum. The shaded area is the 
one-sigma envelope. Sinusoids are either radial velocities 
of the stellar centers of mass (black heavy curve: primary, 
"vradl"; gray heavy curve: secondary, "vrad2"), or of 
points on opposing hemispheres of the primary ("vrotlA" 
and "vrotlB"). Sinusoids are approximate; they do not 
account for the tidal distortions of contact binaries. 

velocity amplitude of a point on the stellar surface is con- 
sistent with the observed perturbation. 

7. Conclusions and Future Work 

The velocity curve of Nex suggests that the emission is 
predominantly associated with the primary (larger, more 
massive) stellar component. This is also consistent with 
coincidence of the amplitude of the count-rate curve with 
the optical amplitude, and with the lack of a distinct pri- 
mary eclipse (smaller, secondary star occulted). The ve- 
locity flip-flop perturbation also suggests that material is 
near the equatorial radius — either compact and near 
the equator itself, or extended at mid- latitudes, but co- 
rotating. 

If the dip in rate seen by IGondoin 20041 in XMM- 
Newton data taken only 10 months earlier is really an 
eclipse, then this indicates rather large and rapid changes 
in coronal structure. The XMM-Newton observation, how- 
ever, did not cover a full continuous period, and the dip 
is offset from phase 0.0, so the interpretation is difficult. 

The line spectrum and crude DEM reconstruction 
show that the the emission comes from a range of tem- 
peratures, and that the emission measure is highly struc- 
tured. The relative abundances of Ne and Fe also appear 
to be similar to other coronal sources: Ne enhanced by a 
factor of a few relative to Fe. 

In the future, we will perform a DEM and abundance 
analysis, in which the line-based simultaneous DEM and 
abundance solutions will be used to predict a better model 
continuum, and thereby improve line flux measurements, 
which are then used to iterate the model. We will also 
perform analysis of composite line profiles in order to re- 
fine the velocity versus phase analysis. By using a dozen 
or so lines, velocity sensitivity can be reduced to about 



30 km s _1 (see |Hoogerwerf, Brickhouse fe Mauche 2004| for 
example). We will perform some analyses of the flare data 
separately to search for line shifts or flare occultation. And 
finally, we will do careful measurement and modeling of 
the density-sensitive He-like triplet ratios. 
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